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Abstract 

Recently, Martin Hutzenthaler pointed out that the explicit Euler method fails to converge 
strongly to the exact solution of a stochastic differential equation (SDE) with superlinearly 
growing and globally one sided Lipschitz drift coefficient. Afterwards, he proposed an explicit 
and easily implementable Euler method, i.e tamed Euler method, for such an SDE and showed 
that this method converges strongly with order of one half. In this paper, we use the tamed 
Euler method to solve the stochastic differential equations with piecewise continuous arguments 
(SEPCAs) with superlinearly growing coefficients and prove that this method is convergent with 
strong order one half. 


Keywords: Stochastic differential equations with piecewise continuous arguments; tamed Euler 
scheme; Strong convergence. 


1 Introduction 

Differential equations with piecewise continuous arguments(EPCAs) represent a hybrid of con¬ 
tinuous and discrete dynamical systems and combine the properties of both differential and dif¬ 
ference equations. They provide a mathematical modeling for physical and biological systems in 
which the rate of the change of these systems depends on its past state. This kind of equations 
plays an important role in many branches of science and industry such as physics, biology and 
control theory, and has been initialed in [24||25| . The general theory and basic results for EPCAs 
have by now been thoroughly investigated in the book of Wiener [26] . 

Since many EPCAs can’t be solved explicitly, computing numerical solutions and analysing 
their properties are necessary. The first work devoted to numerical study for EPCAs is the paper 
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of Liu et al [9]. Subsequently, Song et al [18] and Yang et al m studies the stability of 6- 
methods for advanced equations and Runge-Kutta method for retarded equations, respectively. Liu 
and Gao urn gave the conditions under which the Runge-Kutta methods preserve the oscillation 
of linear EPCAs, and afterwards. Some authors considered the stability and oscillation of the 
numerical solutions [Ml 1221125] , Song and Liu HZ] constructed the convergence improved linear 
multistep method for EPCAs. 

However, systems are often influenced by environmental or some occasional events, which leads 
that the deterministic differential equations can’t demonstrate the real world. In order to avoid 
this problem, many researchers turn to study SDEs. And there have been lots of results on both 
analytical and numerical solutions. The explicit Euler method (see mmm) is most commonly 
used for approximating SDEs with global Lipschitz continuous coefficients. Unfortunately, the 
coefficients of large number of SDEs don’t satisfy global Lipschitz condition Bnamaiaiiaiiaiaa 
.211:12 k 1 . Higham et al [5J showed that the explicit Euler method is convergent strongly when the 
coefficients of SDEs are local Lipschitz continuous and the p moment of both exact and numerical 
solutions are bounded. But, Martin Hutzenthaler in [2] proved that the explicit Euler method 
does not converge in the strong mean square sense to the exact solution of SDE with superlinearly 
growing and globally one-side Lipschitz continuous drift coefficient. To overcome this difficulty, he 
in {3J proposed a modified explicit Euler method, i.e. tamed Euler method in which the drift term 
is modified such that it is uniformly bounded, which is convergent strongly for such SDE. 

Up to now, only a few people considered SEPCAs. Zhang and Song in [29] investigated the 
strong convergence of explicit Euler method for SEPCAs when the coefficients are globally Lipschitz 
continuous or grow at most linearly. Moreover, Song and Zhang in m proved the convergence in 
probability of explicit Euler method under the local Lipschitz and Khasminskii-type conditions. 

Throughout the whole paper, we investigate the numerical solution of the tamed Euler method 
to SEPCAs of the form 

dx{t) = p(x(t), x([t]))dt + a(x(t), x([t]))dB(t), t E [0,T], (1.1) 

with initial value x(0) = £, where B(t) is a r-dimensional Brownian motion, x(t) € M rf , p : —X 

is a continuously differentiable function, a : R d xl^ R dxr satisfies global Lipschitz condition. 
[•] denotes the greatest-integer function. 

2 Notations and main theorem 

Let (D, J 7 , P) be a completed probability space with a filtration {Ft}t >o satisfying the usual condi¬ 
tions and B(t) be a r-dimensional Brownian motion defined on this probability space. Furthermore, 
for any real valued (JZ}" a dapted process x(t), we use ||^(t)||L p (r 2 ,R) to denote (E||x(t)|| p ) p . 

We use notations ||x|| := {x\ + x\ + ... + 2 ^) 2 , (x,y) '■= x± ■ yi + X 2 ■ 1/2 + ■ ■ ■ + %d ' Vd for 
all x = (xi,x 2 ,...,x d ) € R d , y = (yi, y 2 , ■ ■ ■, yd) € R d , and ||A|| := sup a , eK d > || a .||< 1 ||Ax|| for all 
A e M rxd . What’s more, if A is a matrix or vector, its transpose is defined by A T . Set Y17=u = 0 
if u > n. In the whole paper, we make the following assumptions on the SEPCAs (11.11) . 


2 


Assumption 2.1. Assume there exist non-negative constants K,c such that for any x, y, xi, y±, x 2 , 
y 2 € M. d , the coefficients p, and a satisfy 


- cr(a?2,2/2)|| < K(\\xi - x 2 \\ + \\yi - 2/2II), ( 2 . 1 ) 

(xi - x 2 ,p(x 1 ,y) - p(x 2 ,y)) < K\\xi - x 2 \\ 2 , (2.2) 

\\y{x,yi) ~ p(x,y 2 )\\ < I<\\yi - y 2 \\, (2.3) 

\\/J>x(x,y)\\ < K{ 1 + ||a:|| c ). (2.4) 


Remark 2.1. The conditions (12.11) . (12.21) . (12.31) imply that there exists unique solution to equation 
mu>, and the p moments of this solution is bounded in any finite interval [0,T] (see Ifff) j/1. 


In the rest of this paper, let h = ^ be given stepsize with integer m > 1. Grid points t n are 
defined as t n = nh , n = 0,1, • • • . For simplicity, we assume T = Nh, then N = Tm. We consider 
the explicit tamed method for (II.ID . which is defined by taking y 0 = x(0) and, generally 


hidim y\nh\m)h 

Vn+1 VU ^ + h\\n{yn,y[nh]m)\\ 


A @{d)m y[nh\m) 


(2.5) 


where A B n = B(t n+ 1 ) — B(t n ), y n is the approximation to x(t n ). 

In order to formulate the convergence theorem for the tamed Euler method (12.51) . we now 
introduce appropriate time continuous interpolations of the time discrete numerical approximations 
(|2.5I) . More formally, let y(t) : [0, T] x 11 —>• Wf be a sequence of stochastic process given by 


, [* TiVn, y[nh\m) j \ C ( U D/ ^ 

y{t)=y n + I J ,, . - \u ds + I a yyni y[nh]m)dB{s) 


( 2 . 6 ) 


for all t € [t n , t n+ 1 ), n € {0,1,... ,Tm — 1} and m € N = {1, 2,...}. It is easy to get that y(t) is 
{ J-j } — adapted stochastic process. Now we can establish the main theorem of this paper. 


Theorem 2.1. Suppose Assumption \2.1\ hold. Then for each p > 1 there exists non-negative 
constant C dependent on p, K , r and c, but independent of h, such that 


IE 


sup \\x(t) - y(t)\\ p 
te[o,T] 


1 

V 


< Chk. 


(2.7) 


Here x(t) denotes the exact solution of (11.11) . y(t) is the continuous approximation solution of 
tamed-Euler method. 


The detailed proof of this theorem is given in Section 3. 

Remark 2.2. Theorem 1 2. 1\ illustrates that the time continuous tamed-Euler approximations (E2D 
converge in the strong L p -sense with the supremum over the time interval [0,T] inside the expecta¬ 
tion to the exact solution of (ED with the standard convergence order 0.5. 
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3 Estimation of p-moments and proof of Theorem 12.1 

First of all, we introduce several notations here 


and 


A = (1 + T + 2||/r(0,0)|| + 4 K + 211(7(0,0)||) 4 , 


= 1 {||y [ „h] T „||<l l ||Wn||<l}‘ 


Vn 


^il/m y[nh\m) \ n \ 
112 , . 112 / 


^/||2/n|P “1“ 112/^-11 ^ Ill/[n/i] 


n— 1 


D n = (A+ ||£|| + l)exp (4T\/A + 2 sup V' (\/A||ABj|| 2 + a*)), 

«e { -i.o.i.-,n-i} i=u+1 


Tin — \ uj e Tl 


sup Dj (a;) < m ^ , sup 

je{o,i,--- ,n-i} je{o,i,- ,n-i} 


\ABj\\ < 1 


}• 


(3.1) 

(3.2) 

(3.3) 

(3.4) 


here P n € T for all n € {l,2,...,Tm} and m € N. In the following, we give some lemmas which 
will be useful to prove the main theorem 12.11 

Lemma 3.1. If Assumvtion [KT\ holds, and y n ,D n and Tl n are defined by (12.51) . (13.31) . (13.41) . Then 


lf2 n II Vn II 5s D n 


(3.5) 


for all n € (0,1,... ,Tm} and m € N. 

Proof. According to the definition of Tl n , we can see || A5 n || < 1 on f2 n+ i for all n € {0, 1, ... ,Tm— 
1} and mG N. So the assumption (12.11) implies that 


11^71+11| <||2/n|| + h\\n(y n ,y[ nh]m )\\ + \\a(y n ,y [nh]m )AB n \\ 

A || Vn II + hK( 1 + ||2m|| c )||2/?i|| + hK\\y[ nh j m \\ + h||/x(0,0) || 
+ K(\\y n \\ + ||?/[ ri / l ] m ||)||Al? ri || + ||cr(0,0)|| || AB n 11 
<1 + 3 TK + T||/r(0, 0)|| + 2 K + ||cr(0,0)|| < A 


on P n+ i D {||y n || < 1, ||?/[n/i]mll A 1} for all n € (0,1 ,... ,Tm — 1} and m € N. Moreover, the 
estimate 2 ab < a 2 + b 2 for all a, b € M show that 


\\y n+ i \\ 2 =hn \\ 2 + 


\\h(yn,y[nh]m)\\ 2 h 2 


+ 2 


(1 + h\\fj,(y n ,y[ nh]r 
hil/m y[nh]m)h 


+ \W(yn,y[nh]m)AB n 


r j ®{ym y[nh]m)AB r 


i + h\\y(y n ,y[nh]7 

I o/ y{yn, y[nh\m)h / U n 

T2(y n , || , ... + <r(y n , yj n /j] m ) Ai?^ 

<ll2/n|| 2 + 2 \\y(yn,y[nh]m)\\ 2 h 2 + 2||cr(y n , y [n/l]m ) || 2 
(fini h J (Un-. y[nh]nii)tl 


IA B„ 


+ 2 - 


+ 2 ( 2 / 71 , cr(y n , y [nft]m )AS n ) 


1 h||/r(y n , y^ n hj r 

on Tl for all n e { 0 , 1,... , Tm — 1} and all me N. The condition ( 12 . 31 ) . ( 12 . 41 ) implies that 


(3.7) 
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(3.8) 


\\fj,(x,y)\\ 2 <3\\y(x, y) — /i(0, y)\\ 2 + 3||/x(0, y) — /x(0,0)|| 2 + 3||/r(0,0)|| 2 
<3A 2 (1 + ||a;|| c ) 2 ||x|| 2 + 3K 2 \\y\\ 2 + 3||^(0,0)|| 2 
<3AT 2 (1 + m^) 2 ||z|| 2 + 3AT 2 ||y|| 2 + 3\\fx(0, 0)|| 2 
<m(15A' 2 + 3||/i(0,0)|| 2 )(||x|| 2 + \\y\\ 2 ) 

<mV A(||x|| 2 + ||y|| 2 ) 

for all x,y € with {||x|| < 1, ||y|| < 1} C fl {||x|| < m^, ||y|| < m^}. Using the condition (12.11) . 
the following estimate is obtained. 

\\a(x,y)\\ 2 < (|| <r(x,y) — cr(0,0)|| + ||<r(0,0)||) 2 
<(K\\x\\+K\\y\\ + ||a(0,0)||) 2 

< 2(/\ + ||cr(0,0)||) 2 (||x|| 2 + ||y|| 2 ) 

< V\(\\x \\ 2 + ||y|| 2 ) 

for all x,y G R d with {||x|| < 1, ||y|| < 1} C . Moreover, condition (12.21) yields that 
(x,y(x,y)) <A'||z|| 2 + ||s||||/i(0,y)|| 

<A'||x|| 2 + i|M| 2 + A' 2 ||y|| 2 + M0,0)|| 2 
<(K 2 + AT + 1 + ||//(0,0)|| 2 )(||x|| 2 + ||y|| 2 ) 

<\/A(||x|| 2 + ||y|| 2 ) 

for all x,y G with {||x|| < 1, ||y|| < 1} C . Combining (|3.8I) . (13.91) and (13.101) . (|3.7I) yields that 

||2/n+i|| 2 < \\yn\V + (4/i\/A + 2\/X|| AA n ||" + 2a n )(||y n || 2 + ||y[ n /i] m || 2 ) (3-11) 

on {w € 9. : ||y n || < 1, \\y[ nh ] m \\ < 1} C D {w € 9 : ||y n || < ||y[ n fc]mll < for all n G 

{0,1, ■ ■ ■ , Tm — 1} and m € N. 

Before proving (13.51) . we define the mapping 

T n (u) = max{{—1} U{i £ {0,1,2,--- ,n}, ||yi(w)|| < l,\\y [ih]m (u>)\\ < 1}} (3.12) 

for all co G 9, n € {0,1,..., Tm}. With the estimates (13.61) and (13.111) at hand, we now prove (13.51) 
by induction. The base case n = 0 is trivial. Now, let n € {0,1,2, ■■■ ,Tm — 1} be fixed and 
arbitrary. Assume inequality (13.51) hold for all j G {0,1, 2, ■ ■ ■ , n}. Then we show that inequality 
(13.51) holds for j = n + 1, that is 

llj/n+iMH < D n+1 (u) (3.13) 

for all co € Let co € fl n +i be arbitrary, because of the definition of f l n , we can get u ; € 

9 n+ i C f lj which indicates that ||y J (u;)|| < Dj(uo) < m£ for all j G {0,1,2,--- , n}, and it also 
follows from (|3.12l) that 

1 < max{||yj(a;)||, ||y [i/l] (w)||} < m ^ 
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for i = r n (ui) + 1, r n (w) + 2, n. 

Case 1: if r n ( u) > [nh]m, then from (13.111) 

||y n+ i(w)|| 2 + \\y [nh]m (u)\\ 2 

<(l + AhV\ + 2\/A||AS ri || 2 + 2 a n ) (||y n (o;)|| 2 + ||y[ n ft]m( a; )l| 2 ) 

<(ll2/n(w)|| 2 + \\y [nh]m (uj)\\ 2 ) exp (a1iV\ + 2VA|| AB n (uj) || 2 + 2a n (w)) 

< ... 

n 

<(l|yr n +i(w)|| 2 + \\y[nh]m(u)\\ 2 ) exp (4hV\(n - r n ) + 2 (\/A||ASi(a;)|| 2 + <*(<*;))) 

*=Tn + l 

Due to ||y Tn (w)|| < A+ ||f||, || y[ n h]m(u)\\ < 1, we get 

n 

||2/n + iM|| 2 + ||%A]mM|| 2 < (A+ ||£|| + l) 2 exp ^4h\/A(n - r n ) + 2 ^ (VA||ASi(w)|| 2 + ^(w))). 

i=r n +l 

(3.14) 


Therefore 

l|j/n+l(w)|| < 



sup 

u£{—1,0,1 


.,n} 


(V\\\ABi(uj)\\ 2 + Q! i (w))^ <D n+1 (uj). 

i=u -\-1 


(3.15) 


Case 2: if t„(w) < [nh]m, Using 2ab < a 2 + 6 2 , we derive ||y n +i|| from (13.111) and (13.141) . and 
obtain 


^Hl/n+l (k-0 || X ||y[nfe]m(^) II — H^ra-fl (^) || T ||y[n/i]m(^) II "" 

n 

<2||y[ n / J ] m (w)|| 2 exp (AhV\(n - [nh\m) + 2 ^ (\/A||ABj(w)|| 2 + a;(w))). 

i=[nh]m 


Hence 

n 

Hs/n+iMII < \\y[nh]m(u)\\ exp (Wa/i(ti - [ nhjm ) + 2 ^2 (VA||AHj(w)|| 2 + a*(u;))) 

i=[nh]m 


(3.16) 


< ... 


< ||y[r„h+i]mM|| exp (4hV\(n - [r n h + 1 \m) + 2 ^ (VA||AHj(u;)|| 2 + «*(<*;))) 

i=[T n h+l]m 

n 

< (lbr„+i(w)|| + ||y[ Tn ft] m ||) exp (4hV\(n - r n ) + 2 XZ {V\\\ABi(uj)\\ 2 + ai(u)) 

i=[T„h+l}m 


< 


[T„h+l]m 

+ XZ (v / A||A.B i (u;)|| 2 + ai(w))) 

i=T n + l 


+ 2 sup 

«€{ —1,0,1, ...,n} 


X] (VAIIA^Mf + a^u;)))) <D n+1 ( w ). 

2=li+l 


(3.17) 
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The proof is completed. □ 

The following two lemmas are useful to prove that D n is bounded on 17. 

Lemma 3.2. J3j/ For all p > 1, 


Tm—l 


sup E 

m>4Ap 


exp (p\ 22 ll A ^H 2 ) 


< e 


2 A pTr 


(3.18) 


«=o 


Lemma 3.3. Let on : 17 —>• M for all i € {0,1,... ,Tm}, and m € N be given by (13.21) . Then for 
each p > 1, 

n —1 

l|2\ 


sup sup E 
zG{—1,1} m£N 


sup exp^Vai) < 2e 2pT ^ + ll <T (0’0)ll ). 
nG{0,l,...,Tm} 


(3.19) 


i =0 


Proof. Note that the time discrete stochastic process z ^"Tq 1 is an artingale for every 

n € {0,1,... , Tm}, z € {—1,1} and every m € N. Then it is easy to deduce that the time discrete 
stochastic process exp (z Q *) is a positive Tt n -submartingale for every n € {0,1,..., Tm}, 

z € {—1,1} and every m € N. Hence Doob’s martingale inequality shows that 


71—1 

sup exp (z af) 
ne{0,l,...,Tm} i=0 


< 


P 


L p (Q,R) p — 1 


Tm—l 


exp (z 22 “0 


i =0 


Z/p(Q,]R.) 


(3.20) 


for all m £ N, p G (1, +oo) and all z € {—1, !}• Moreover, we have that 


E 


M{IMI<i,NI<i} c (- 


°(x,y) 


A 5i)| : 


\/IMI 2 + IMI 2 ’ vlMF-HM 2 

/ 2 ui \\x T °{x-,y)\\ 2 

<P h 1 { N |2<i iW2 <i}c^^||2 + || y ||2)2 

2(A' + ||<r(0, 0)|| 2 )||x|| 2 

<P /ll{||*||a<l,||v||a<l}c-n^p + || y || 2 - 

<2p 2 h(K + || cr(0,0) || ) 2 

for all x € M rf , i € {0,1, ...,Tm — 1}, m € N, p £ [l,oo) and all z £ {—1,1}. Therefore lemma 4.3 
in [1J gives that 

x cr(x, y) 


E 


exp (p~l{||x||<i,||y||<i} c ( 


V \\ x \\ 2 + IMI 2 ’ OTF 


ASi)) 


(3.21) 


< exp (2 p 2 h(K + ||cr(0,0)||) 2 ) 


for all x € M d , i € {0,l,...,Tm — 1}, m £ N, p € [l,oo) and all z € {—1,1}. In particular, (13.21[) 
shows that 

E exp(pzdj)|< exp (2p 2 h(K + ||ct( 0, 0)||) 2 ) 
for all i € {0,1, ...,Tm — 1}, m € N, p € [l,oo) and all z € (—1,1}. Hence, we obtain that 

Tm—l 

E exp (jpz E a,)] < {f2p 2 T(K + ||<r(0,0)||) 2 ) (3.22) 

2=0 
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for all m € N, p £ [1, oo) and all z £ {—1,1}. Combining ()3.20f) and (|3.22l) . then for all p £ [2, oo) 

< 2 exp (2 pT(K + ||<r(0,0)||) 2 ). 


sup sup 

ze{-i,i}meH 


n— 1 

sup exp (z a i ) 
ne{0,l,...,Tm} i=Q 


L p (n, R) 


The proof is completed. 

Lemma 3.4. Let D n is defined by (I -l.3 1) . Then for all p £ [l,oo) 


sup E 


sup \D n \ p 


< oo. 


mG N l 0 <n<Tm 

Proof. The proof is similar to that of lemma 3.5 in [3j. 

Lemma 3.5. Let LtTm £ J~ for m £ N be given by (|3.4jh Then for each p > 1 we have 

sup (m p ■ P[(n Tm ) c ]) < oo. 
raG N 


□ 

(3.23) 

□ 

(3.24) 


Proof. This Lemma is based on the Lemma 13.31 and 13.41 and the proof is same as that of Lemma 
3.6 in [3]. □ 

Next, we will prove the boundedness of y n in L p sense. 

Theorem 3.6. Let y n : —>• M d , for n € {0,1, • • • , Tm} and m € N be given by (|2.5p . then for all 

p G [1, oo) 

1 < oo. (3.25) 


sup 

meM 


sup E||y n || p 

0 <n<Tm 


Proof. First, we can by (12.511 represent the approximation y n as following 

n —1 r / \ n —1 

Vn =£ + cr(0, 0 )B n + ^2 i jn J ) V[lh]m + 22 V[ih]m ) ~ ^(0, 0)) AS* 

for all n € {0,1,... ,Tm}, and m € N. The Lemma 4.7 in [T] and Burkholder-Davis-Gundy type 
inequality in Lemma 3.8 in [3] then give that 

II2/nlli/p(o,R) <||f||z, p (fi,R) + pVrT\\cr( 0 , 0)|| +Tm 

n —1 r i 

+ hv(yiiy[ih]m) — cr^(°, °)|li, 1J (n j R)/i 

i =0 v=\ 

<(HCl|Lp(n,R) + pVrT\\cr(0, 0)|| + Tm ) 

n —1 


+ pKV5rh(22 (lli«lli p (n I R) + \\y[ih]m\\l p (n,i 

i=0 


Using Gronwall’s inequality, we can get 

llynlli p (n ) R) < 2 (ll?llLp(n,R) +pv'rf||<r(0,0)|| +Tm)‘ 


n— 1 


Sill p (n,R) 


+ 8 p 2 K 2 rh22 SU P 11%' 

i=o 

< 2 (IICIlLp(a,R) +pVPT\\(T(0,0)\\ +Tm) 2 e 8p2K2Tr 
















for all n G {0,1,... ,Tm},m G N and p G [2, oo). For all n G {0,1,2, m G N and p G [2, oo), 

it is easy to obtain that 


l|yn||L p( n,R) < >/2(||£|| M n,R) +pV?T\\a(0,0)\\+Tm)e 4 P 2K2Tr . (3.26) 

Of course (13.261) doesn’t prove ||yn||L p (n,R) < oo, due to m G N on the right-hand side of (13.261) . 
However, Holder inequality and lemma 3.5 show that 

sup sup ||l(rc„)c2M||z, p (n,R) 
mgN nG{0,l,...,Tm} 


(3.27) 


< sup sup (m||l(n„)c||L 2 p(fi,]R)) x ( sup sup m 1 \\y n \\ L2p (n,] 

m £N ng{0,l,...,Tm} ' mgN ng{0,l,...,Tm} 

<V2e 4p2 ^ 2Tr ( sup m 2p • P[(H Tm ) c ]) 2p x (||^||z, 2p (r2,K) + pVrT\\a(0, 0) || + T) 

'' mgN ' 


<oo 


for all p G [2,oo). Additionally, Lemmas 13. II and 13.41 give that 

sup sup ||ln„yn|U p (n,R) < sup sup ||£> n || Lp(n>R) < oo 
mgl n£{0,l,...,Tm} mgN n£{0,l,.,.,Tm} 

for all p G [2,oo). Combining (13.271) and (13.281) the theorem is proved. 


(3.28) 


□ 


Lemma 3.7. Let y n : Q —>• U. d , for uGN = {1,2, ..., Tm} and m G N be given by (|2.5I) . then for 
all p G [1, oo) 

sup sup m\p(yn,y[nh]m)\\ P ] < oo, sup sup E[||er(y n , y [nh]m ) || p ] < oo. (3.29) 

mgNn£{0,l,...,Tm} mgM ng{0,l,...,Tm} ' 


Proof. By condition ()2.41) . we can get for any x,y G 


\\h{x,y)\\ <11^(0, 0)|| + ||/z a: (0i®,y)||||®|| +K\\y\\ 
<11^(0,0)11 +K( 1 + ||a;|| c )||a;|| + K\\y\\ 
<|| M (0,0)|| +2^1 + ||z|| c+1 ) + iqy||. 


It comes from theorem 13.61 that 


sup sup \\p(y n , y[nh]m) llz, p (fi,R) < 11^(0, 0)|| + 2AT I 1 + SUp SUp 
mgN ng{0,l,...,Tm} ' m£Nng{0,l,...,Tm} 


U 


e+1 \ 


A K SUp SUp ||y[nh]mlU p (n,R) < OO 

mgN n£{0,l,...,Tm} 


for all p G [l,oo). According to the global-Lipschiz condition (12.11) and theorem 13.61 for any 
x, y G K rf , we obtain 


sup sup \W(yn,y[nh]m)\\L p (n,R) < |]cr(0, 0) || +K ( SUp SUp | |l/r» 11 («,R) ) 

m£Nn£{0,l,...,rm} 'm£Nn£{0,l,...,rm} ' 

+ AH SUp SUp ||y[n/i]mlll/ p (S7,R) ) < OO. 

'm£Nii£{0,l,...,Tm} ' 

□ 
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Now we are in a position to give the proof of theorem 12.11 


Proof. We define t = t n for any t G [t n , t n+ \), and n G {0,1,... , Tm — 1}, m G N. It is known from 

dMD 

rt 


y(t) = £ + f 

Jo 

for all t G [0, T]. Note that 


o l + h||/r(?/(s),y([s]))|| J 0 

V(y(s),y([s})) 


ds+ [ cr(y(s),y(\s]))dB(s) (3.30) 

Jo 



l + h\\y(y(s),y([s]))\\ 


ds 


for all t G [0, T] P-a.s.. Hence, Ito’s formula yields that 

\\x(t) - y(t)\\ 2 = 2 [ ( x ( s ) - y(s),n(x(s),x([s])) - y{y(s),y{\s])))ds 
Jo 


+ 2h f („ w - vW , 

+ I ||o'i/(*(s),*([s]))-o'«/(y(s).2/([«]))|| 2ds 

v=l J ° 

J rt 

+ 2 


!E / <x(s) - ?/(s),cr [/ (x(s),x([s])) - cr I/ (y(s),y([s]))>dH ;y (s) 
„=i - y o 


=Ai + H 2 + H 3 + H 4 


for Hi, H 2 , H 3 , we can estimate 

Hi =2 / (x(s) - y(s), it(x(s), x([s])) - y(y(s), x([s])))ds 
Jo 

+ 2 / <x(s)-y(s),//(y(s),x([s]))-//(y(s),y([s])))ds 
Jo 

+ 2 [ (x(s) - y(s),n(y{s),y([s])) - y(y(s),y([s])))ds (3.31) 

Jo 

<(2K + 2)f \\x(s) - y(s)\\ 2 ds + K 2 [ ||x([s]) - y([s])|| 2 ds 
jo Jo 

+ [ ||My(s),y([s])) - y{y(s),y{[s]))\\ 2 ds, 

Jo 

A 2 < [ \\x(s)-y{s)\\ 2 ds + h 2 f \\y{y(s), y([s]))|| 4 ds, (3.32) 

jo Jo 

H 3 <2V f ||av(x(s),x([s]) - cr u (y(s),y(\s}))\\ 2 ds 
u=i J o 

+ 2 Y1 [ || cr ^(?/(s),2/([s]))) -a u {y(s),y([s]))\\ 2 ds (3-33) 

u=l J 0 

rt rt 

<4 K 2 r / \\x(s) — y(s)\\ 2 ds + 4K 2 r / ||a;([s]) — y([s])|| 2 ds 
Jo Jo 
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+2 K 2 r 


[ ||l/(s) — 2 /(s)|| 2 o?s- 
J 0 


The Global-inequality, (13.31 p , (13.321) , (13.33P and 2 ab < a 2 + b 2 shows that 

rh 
i o 


sup ||x(f) — 
te[o,ti] 


m 

< {AK 2 r + 2K + 3) / ||x(s) — y(s)\\ 2 ds 

Jo 

+ (4r + 1)K 2 f \\x([s})-y([s})\\ 2 ds + h 2 [ \\n{y(s), y([s]))|| 4 ds 
Jo Jo 


+ 


[ ||/4y(s),Z/([s])) - v(y(s),y(\s}))\\ 2 ds + 2K 2 r f \\y(s)-y{s)\\ 2 ds 
Jo Jo 

r rt 


(3.34) 


*- is./ 

tE[0,ti] i/=\ ^0 


+ 2 sup 


(x(s) - y(s), a v (x(s), x([s])) - a u {y(s), y{\s])))dB u (s) 


P-a.s. for all t\ G [0, T]. The Burkholder-Davis-Gundy type inequality in Lemma 3.7 in [3j hence 
yields that 


sup \\x{t) - y{t)\[ 
te[o,ti] 


ip(nR) < J o ||Mi/(a),I/([s]))-My(s)»I/([s]))|L p (Q,M) 


I Lp (Q,]R) 


+ (8K 2 r + K 2 + 2K + 3) f ( sup || x(u) — y(u) 

Jo ^ue[o,s] 

f T f T (3.35) 

+ 2iL 2 r \\y{s) - y(s)\\ 2 Lp{nM) ds + h 2 \\y(y(s), 2/([s]))||l 2p(n , R) ds 

+ 2 p (S/“ ||(*(s) - 2 /(s),a„(x(s),x([s]) -a v (y( 5 ) 1 y(fe])))||l J , (niR) *) a 


for all t\ G [0, T], and all p G [4, oo). Next the Cauchy-Schwarz inequality, Holder inequality and 
again the inequality 2 ab < a 2 + b 2 for all a, b G M imply that 


a 


^ 2 p{Y1 [ |*(a)-y(a)||i p (n 1 R)lk«'(*(»)»*([a])) -^(l/U)»«(|s]))|li 1I (n I R) ds j 

u=l J0 

sup ||x(s) - y{s) || M n )R) ) (X) [ ||<r„(x(s),x([a])) - ^( 2 /(s),y([s]))||^ (n]R) ds) 2 

sE[0,£iJ ,,_-| J 0 


rt\ 

— 2 s r u P , IK*) " 2/WllL(n,M) + A P 2R2r / (IM S ) “ yU)lli P (n,R) + IHW) “ v(\s\ )lli p (n,R)) ds 


te[o,ti] 


< 


1 

2 tg s JP IK*) “ + 8 P 2R2r l IK S ) - y( s )lli p (n,R) ds 

+ 8 P 2 K 2 rj^ \\y(s)-y(s)\\ 2 Lpm ds + Ap 2 K 2 r ^ ||x([s]) - 2/( W)|| 


(3.36) 


for all t G [0, T], and all p G [4, oo). Inserting inequality above into (13.351) and applying the estimate 
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(a + b ) 2 < 2a 2 + 2 b 2 for all a, b G M, then yield that 

sup ||x(t)-y(t)|||| <(2K 2 r + 8p 2 K 2 r) f ||y(s) - j/(s)||* p( n,R) ds 

+ (8K 2 r + K 2 + 2K + 3 + \2p 2 K 2 r) f ( sup || x(u) — y(u 

Jo ' ugfO.sl 


I L/p (f2,K) 


(3.37) 


+ k2 L H^-)’ y ([-]))ll£2 P m ds + J 0 llM2/(s),2/([s)) -/ x (yU)>y(fe]))llLp(n,]R) ds 

for all t G [0, T], and all p G [4, oo). In the next step Gronwall’s Lemma shows that 


sup \\x(t) - y(t)\\ 
te[o,ii] 


Lp 

T 


< 2 e T{ - SK2r+ K 2 + 2K +z+Wp 2 K 2 r) x (^2K 2 r + 8p 2 K 2 r) 


+ 


J o IbW - »(l)|lh (a i i) da + h 2 f 0 W»(s).»(W))|| 
f \\Mv(s),s(m- p (Q,M) 

J 0 


L2 P (r2,R) 


ds (3.38) 


and hence, the inequality y/a + b + c < y/a + y/b + y/c for all a, b, c G [0, oo) gives that 


sup ||x(t) — 
te[o,ti] 


Lp (Q,]R) 


<V 2 T e h T ( 8K2r +K 2 +^ K +^p 2K2 r) x ( R 2 r + 8 p 2 K 2 r 

+ h 


sup 

te[o,T\ 


Lp (r2,M) 


(3.39) 


sup \\p(y(t),y{[t}))\\ ( ) + sup \\p(y(t),y(\t]))-p(y(t),y(\t}))\\ , } 

L te[o,T] pK ’ n te[o,T] py ’ ’ 


for all N G N and all p G [4, oo). Additionally, the Burkholder-Davis-Gundy inequality in Lemma 
3.7 in P| shows that 


Lp (Q,R) 


sup ||y(t) -y(i)|| Lp(f2iR) < sup 


te[o,T] 


t£[0,T] 


+ sup || f u(y(s),y([s]))dB(s)|| , R) 

te[o,T] Jt pv ’ ) 


<h( SUp \\p(y n ,y[nh}r 
' n£{0,l,...,Tm} 


L/p 


t r 


+ P 


r u =i 


<h( SUp || d{yn,y[nh]m 

' nG{0,l 


Lp (Q,M) 


+ pVhr( SUp SUp || <7v(yn,y[nh\r 

' z/G{0,l,...,r} 7iG{0,l,...,Tra} 


L p (Ll,M) 


for all m G N and all p G [2,oo). Lemma [3771 implies that, for all p G [l,oo), there exists C\ > 0 
independent of h such that 


sup 

te[o,T] 


Lp 


< Cih5. 


(3.40) 
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In particular, we obtain that for all p £ [l,oo) 


Moreover, the estimate 


sup ||2/(£)|| L (QK) 
iG[0,T] PV ’ ; 


< OO. 


(3.41) 


sup \\p(y(t),y([t})) - p{y(t),y([t]))\\ fnm 
te[o,T] v\ , ) 



L/2cp{Q^) 


sup \\y(t) - y(t)\\ 

ie[0,T] 


L2p(^2,R) 


(3.42) 


for all p £ [1, oo). Then inequality (13.401) and (13.411) hence show that there exists C *2 > 0 independent 
of h 


sup \\p(y(t),y([t])) - n(y(t),y([t}))\\ L (nm <C 2 h* 
te[o,T] M 1 


(3.43) 


for allp € [1, oo). We obtain from (I3.39p . (|3.40l) . (13.43|) and lemma lTTl that there exists non-negative 
constant C independent of h such that (12.71) holds. The proof is complete. □ 


4 Numerical Experiments 

In this section, we give two numerical experiments to illustrate the strong convergence and the 
convergence order, 
we consider 

f dx(t) = ( — x(t) a + a(x(t) + x([t])))dt + b(x{t ) + x([t]))dB(t), t £ [0,T], ^ 

|^x(0) = c. 

In the first numerical experiment we used the parameters a = 3, a = 0.5, b = 1, c = 1.5. In 
the second numerical experiment, we used parameters a = 5, a = 4.5, b = 3, c = 1. 

We square both sides of (12.71) with p = 2, we get the mean square error E[sup ig r 0 T i || x(t) — 
y(t)|| 2 ] which should be bounded by Ch. The mean square error at time T was estimated in the 
following way. A set of 30 blocks each containing 100 outcomes(w^ : 1 < i < 30, 1 < j < 100) were 
simulated. We denoted by y(T,uoij ) the numerical solution of the jth trajectory in the ith blocks 
and x(T,ujij ) the exact solution of (14.11) in the jth trajectory and ith block. The ’exact solution’ 
was computed on a very fine mesh(we used 262144 step). 

Let e denote the mean square error. Then by the law of large numbers, we conclude that 

1 30 100 

< T ) = ^EE - y(T,uij Oil 2 - 

*= 1 3 =1 

There are three test in each numerical experiment with T = 1,2,3. We can see from the table 1 
and table 2, the ratios of errors in the tables are consistent with the theoretical rate of convergence 
as stated in theorem O 
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step 

e(l) 

ratio 

e(2) 

ratio 

6(3) 

ratio 

2-8 

0.0022 

* 

0.0038 

* 

0.0089 

* 

2—9 

0.0010 

2.2000 

0.0020 

1.9000 

0.0051 

1.7451 

2-iO 

0.0005 

2.0000 

0.0010 

2.0000 

0.0018 

2.8883 

2— 11 

0.0002 

2.5000 

0.0006 

1.6667 

0.0009 

2.0000 

2-!2 

0.0001 

2.0000 

0.0003 

2.0000 

0.0004 

2.2500 


Table 1: The error at times T = 1,2,3, for the first numerical experiment. 


step 

e(l) 

ratio 

e(2) 

ratio 

6(3) 

ratio 

2-8 

0.0379 

* 

0.1550 

* 

0.2779 

* 

2-9 

0.0150 

2.5252 

0.0844 

1.8359 

0.1311 

2.1198 

2— 10 

0.0073 

2.0592 

0.0443 

1.9068 

0.0808 

1.6225 

2— 11 

0.0033 

2.2375 

0.0158 

2.8107 

0.0471 

1.7155 

2— 12 

0.0015 

2.1252 

0.0107 

1.4727 

0.0401 

1.1746 


Table 2: The error at times T = 1, 2, 3, for the second numerical experiment. 
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